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Abstract. To analyze the limiting spectral distribution of some random block- 
matrices, Girko [7] uses a system of canonical equations from [||. In this paper, 
we use the method of moments to give an integral form for the almost sure 
limiting spectral distribution of such matrices. 

1. Introduction and main result 

A random block-matrix is a matrix whose entries are random matrices. In [6], 
Girko studied the spectra of large dimensional random block-matrices by intro- 
ducing a system of equations, called the system of canonical equations, to analyze 
the spectra. This system of canonical equations was used later by Girko [7] to 
study a model for which the system is solvable. The model studied there has many 
restrictive conditions. 

In the current paper, we are going to study the same model under different 
conditions for the blocks. The main tool of the proof is the method of moments. 
We will follow the proof of the main theorem by propositions, as applications to the 
theorem, in which the blocks are made of some known ensembles like the Gaussian 
unitary ensemble and the Wishart random matrix. Free probability theory is used 
to prove these propositions. 

The spectral measure of an n x n Hermitian matrix A is 

1 " 
A*a = - V £a, , 
n 

i=l 

where Ai < A2 < • • • < A„ are the eigenvalues of A. In this paper we consider 
random matrices, i.e., matrices where the entries are random variables on some 
probability space. In this case, /za is, of course, a random measure. 

We will denote the weak convergence of probability measures by lim n — >00 /i n — /1 

or 

v 

fJ-n — ► /i asn-^ 00. 
If the moments of measures converge, 



lim / x k /i n (dx) = / x k [i(dx) 

for all k > 1, we will write 

m 

A*n — > M as n ^ 00. 
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If {/*„} is a sequence of random measures which converges in one of the above senses 
almost surely, we will append the abbreviation "a.s." to the above notation. Wc 
note that for k > 1, the k th moment of /ia is 

(1) J x h ti K {dx) = tr„(A fe ), 

where tr„(A) := i Yh=i A U- 

The Kronecker product ® of two matrices A = (ffly)*j- = i and B = (i>y)*- =1 is 
defined to be the nfc x nfc matrix given by 

oilB ai 2 B ... aifeB 



A <g B = 



a 2 iB a 2 2B . . . a2fcB 



afciB a^B . . . afefcB 

Among the properties of the Kronecker product, we will need the identity 

tr fe „(A <g> B) = tr fc (A)tr n (B). 

Also, if A and C are two k x k matrices and B and D are two nx n matrices, then 
(A <g B)(C <g D) = AC <g BD. Finally, I k is the k x k identity matrix. 
Now we are ready to state the main theorem. 

Theorem 1. For n,k > 1, consider the double array of random block-matrices 
{Bn.fc} whose terms are given by = Ife ® A„ + Wfe (g B„, where for n > 1, the 
matrices A„, B„, and W„ are Hermitian random matrices of order n, and satisfy 
the following hypotheses: 

(i) There exists a compactly supported probability measure fi w such that 

(2) Mw„ —> as n — ► oo a.s. 

(ii) For real t, there exist probability measures ip(t] .) such that 

MA„+tB n — ► ip(t; •) as rt — > oo a.s. 

anrf ^(t; .) has a support that is uniformly bounded for t in any compact 
subset o/K. 

Under these conditions we have 

(3) lim lim [iB„ k = lim lim ^b„ fc = v a.s., 

n — *oo k — >oo ' k — >oo n — >oo ' 

where the probability measure v is defined as 

(4) v(dx) = / ip(t;dx) fi u (dt). 



Remark 1. Since the support of ip(t; dx) is uniformly bounded for t in supp(^) (the 
support of Hu), then the probability measure v(dx), introduced in is compactly 
supported. 

As mentioned above, matrices of the form (g A„ + Wfe (g B n were analyzed 
by Girko [7]- In [3 Theorem 3], Girko assumes that {A„} and {B„} are two 
sequences of real symmetric non-random matrices, B„ is a positive definite matrix 
for each n > 1, and the entries of the symmetric matrix Wfe are independent ±1 
with probability i. He shows under these assumptions that the spectral probability 
distribution Fa n k (x) := fim n k ((— oo, x]) of the sequence of random block-matrices 



LIMITING SPECTRA OF GIRKO'S BLOCK-MATRIX 



3 



{B>n.fc} converges, for almost all x's and with probability one, as both k and n go 
to infinity to a non-random distribution function that follows from a complicated 
equation given in [7] . 

In Theorem [U our assumptions allow us to identify the limit. In the course of 
our proof, we are also able to derive Girko's SS-Law (a sum of semi-circular law), 
see Proposition [H 

2. Proof of Theorem [T] 

We need the following lemma. 

Lemma 1. (5j p. 94, Lemma 3.1] Fix k G N, let T = {to, t\, . . . , t/.} be a set of 
distinct points in R and P n (t) — ao,„ + a\, n t + • ■ ■ + a,k : nt k be a polynomial with 
a^ n G C for every i and n. If P n (t) converges for every t G T as n — > oo, then 
the limit is a polynomial of degree < k, say it is P(t) — oq + ait + • • • + akt k . 
Moreover, the convergence is uniform on every compact subset ofR. Furthermore, 
lim^oo a i<n = a; for every i. 

Proof of Theorem |7J The proof is based on the method of moments. Using 
the aforementioned properties of the Kronecker product, the m th moment of the 
spectral measure of B„.fc is given by 

tr„ fc (B™ fc ) = tr fcn ((I fc ® A n +W k ®B n ) m ) 

(5) 

= E"lo tr fe( W fc) tr„0(A„,B„;m - 

where <f>(A. n , B n ; m — j, j) is the sum of all the noncommutative monomials in which 
B„ appears j times and A„ appears m — j times. Let the j th moment of fj, u be ojj, 
j > 1. By (dJ and (2}, lirm^oo tr fc (W£) = uj a.s. Therefore, 



(6) lim tr rlfe (B™ fc ) = Vw, tr„(0(A„, B„; m - 

k^oo L — ' 
3=0 

On another hand, for all t G M 



(7) tr„((A„+tB„) m ) = £Vtr n (0(A n ,B„;m- 

3=0 

Therefore, 

~ m 

/ tr n ((A n + tB n ) m )/i u ((ft) = y~\uj tr n (0(A n ,B n ;m- 
Jk j= o 

So by ©, 

(8) lim tr nfe (B™ fc ) - / tr„((A„ + tB n ) m )^{dt). 
But, 



lim tr„((A n +tB n ) m ) = f x m ^(t;dx) 

n — >oo Jj^ 



and by Lemma[T]this limit is uniform in t as t varies over the compact set supp(^ 
Therefore, 



lim lim tr n fc( 

n — >oo k — >oo 



U) = J ^x m iP(t;dx) (i u (dt). 
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By Fubini's Theorem 

lim lim tr nfe (B™ fc ) = / x m ( [ i/)(t;dx) fJ, u (dt) 

n->oofc->oo ' J R \J R 

The other iterated limit follows from the observation that ([5]) and j7]) imply 
tr„ fc (B™ fc )= f tr n ((A n + *B n )"> Wfc (di) 

JR 

for every n and k. Since /iw fc is a discrete measure, 
(9) Urn tr nfc (B£ fc ) = / [ x m ^(t;dx)fi Wk (dt). 



n — >oo 



By Lemma U J m x m i^(t;dx) is a polynomial in t and since /iw fc converges in mo- 
ments, it follows that 



lim lim tr Ilfe (B™ fe )= / x m / xp(t;dx) fi u (dt) 

Now, since v{dx) := J R ip(t;dx) ^(dt) has a bounded support, the result follows. 

□ 

Remark 2. The integral in |4} always exists because ip(t; .) is measurable in t. This 
can be seen as follows. The characteristic function of ip(t; .) is analytic, as i/j(t; .) 
has a compact support for each t. So the characteristic function is measurable in t 
as a pointwise limit of the series in the moments f R x k ip(t; dx), k > 1; the latter are 
polynomials in t by Lemma [TJ Therefore the inversion formula of the characteristic 
function implies the measurability of ip(t; (— oo,a]) for any a. 

3. Applications 

In this section, we apply Theorem [T] to some well-studied ensembles of random 
matrices. To do so, we introduce these ensembles and review the pertinent topics 
from free probability theory. 

3.1. Random Matrix Theory. We call annxn Hcrmitian matrix A = (Aij)" J=1 
a Wigner matrix if it is a random matrix whose upper-diagonal entries are indepen- 
dent and identically distributed complex random variables such that E(Aij) = 
and E{\Aij\ 2 ) — — for all i < j. Moreover, the diagonal entries are independent and 
identically distributed real random variables such that E(Au) — and E(A^ i ) = —. 
We will denote all such Wigner matrices of order n by Wigner (n). 

An important example of a Wigner matrix is the Gaussian Wigner matrix for 
which {SL4,j : 1 < i < j < n}U {'SAij : 1 < i < j < n} is a family of independent 
random Gaussian variables such that An ~ N(0, —) for every i and SiA^, 3Aj ~ 
N(0, A-) for every i < j. We will denote all such Gaussian matrices of order n by 
G(n). 

We call the random matrix B = X*X a Wishart matrix if X is a p n X n matrix 
whose entries are complex independent Gaussian random variables such that $lXij , 
QXij ~ N(0, for every 1 < i,j < n. Here X* is the conjugate transpose of X. 
We will denote all such Wishart matrices of order n and shape parameter p n by 
Wishart (n,p n ). See [2111] for more details and references. 
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For these type of random matrices, the limiting spectral distributions are known. 
If A„ is Wigner(ri), then by Wigner's Theorem (c/. [2]), 

v 

/LtA„ — ► 7o,i as 11 ^ oo a.s. 

where ^ a ,a 2 is the semicircular law centered at a and of variance a 1 which is given 
as ^ 

7 Q . (T 2(dx) = j \/ 4cr2 ~ ( X ~ a ) 2 1 [a-2a, a +2a]{x)dx. 

Zncr 

On the other hand, if B„ is Wishart(n,p n ) for all n and lim„^oo — = a > 0, then 
('/• ®) 

v 

/iB„ — * Pa as n — ► oo a.s. 
where p a is the Marchenko-Pastur law with mean a > which is given as 

p a {dx) = (1 — a) + 5o(dx) 

The following two propositions are consequences of Theorem [1] 

Proposition 1. Let {W„}, {A n } and {B„} be three independent sequences of 
random matrices such that W„ is Wigner(n) and A„,B„ are G(n) for all n. 
Then the almost sure limiting spectral distribution v o/B„.fe = Ij, <g) A„ + Wj, ® B„, 
see ([3]), is absolutely continuous with the probability density function 

gi(x) ; whenever 2 < |x| < 2-\/5 



where 



and 



52(2:) ; whenever \x\ < 2 

2 v / 4(l _ n2)— ^2 ^/J— £2 



v^TT (i + t 2 ) 

2 V4(l + i 2 )-a: 2 yi^F 



c/f. 




FIGURE 1 . The probability density function corresponding to the 
limiting spectral distribution of B^fc when A n and B„ are Gauss- 
ian matrices. 

In order to state the following proposition we first define the following functions: 
hi(x;t) = 2 + 27t 2 -3tx -3t 2 x 2 + 2t 3 x 3 , 
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hi (x; t) = 1 — t x + t 2 x 2 , 



H(x;t) = -j= (hi(x;t) + s/hjx; t) 2 - 4 h 2 (x; i) 3 




the two functions si(i) and s^it) which are the two real roots of the quartic equation 
in x 



Proposition 2. Let {W„}, {A„} and {B n } be three independent sequences of ran- 
dom matrices such that W„ is Wigner(n), A„ is G(n) and B„ is Wishart(n,j>„) 
for all n and linin^oo 2a. = 1. Then the almost sure limiting spectral distribution v 
ofV> n ,k = Ife ® A„ + Wfe ® B„, see ([3]), is absolutely continuous with the probability 
density function 



where f(x;t) is given in Equation (| 1 1 [) . 

We simply need to verify the hypothesis of Theorem [TJ We use the machinery 
of free probability (specifically free additive convolution) to do this. 

3.2. Free Additive Convolution. Let /i be a probability measure with a com- 
pact support in R. We define its corresponding Cauchy (Stieltjes) transform to be 
Gn(z) := J R —^fi(dx), for z € C such that ^(z) > 0. The Cauchy transform G^(z) 
possesses the following properties: 

(i) 3(G M (z)) < whenever 3(z) > 0. 

(ii) lim^i^oo zG M (z) = 1. 

(iii) Gfi(z) is analytic in a neighborhood of oo. 
The Stieltjes inversion formula is given by 



for every continuity set E C 0(R) (the a— field of Borel subsets of R). The R- 
transform of /i is defined as R^{z) = K^{z) — - where K^{z) is the inverse function 
of the Cauchy transform G^(z), i.e., G^(K^(z)) = z. The two functions K^(z) and 
Rfj,{z) are well defined in < \z\ < r and < \z\ < r, respectively, for some r > 0. 

The free additive convolution of probability measures with compact supports 
in R arises in free probability theory (c/. |10j). If /i and v are two probability 
measures with compact supports in R, then their free additive convolution EB v is 
a probability measure with a compact support in R, see |10j . The R-transform of 
fiBvis given by R^{z) = R^(z) + R v {z). 

Denote the dilation D t of a measure /x by D t (n), where D t (n)(E) = fi(E/t) for 
every E if t ^ and D t (ji) = 8q if t = 0. Since Ro t {p) i z ) — tR^(tz) for every t e R 
(cf. HOI P-26]), therefore 



(10) 4 + 27 i 2 - 6 t x - x 2 - 6 1 2 x 2 + 2 i x 3 + 4 i 3 x 3 - t 2 x 4 = 

(see Proposition [3] for details) and the probability density function 





(12) 




DtOiBv) = Dt{n)mD t {y). 
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For the limit laws described above, i? 7o ^ 2 (z) — a 2 z and R Pa (z) = (c/. |10j); 
furthermore 70,1 EB D((7o,i) = 7o,i+t 2 for every real t. The next proposition com- 
putes 70,1 EB D t (pi) for £ > 0. 

Proposition 3. If p — 70,1 and v — p\ then for t > 

juffl D t (v){dx) = f(x;t)dx 

where f(x;t) is given by Equation (111)) . Furthermore, for each compact set Ccl 
there exists M > such that the support of f(x;t)dx is contained in [—M,M] for 
all t in C . 

Proof. Fix t > 0. The R-transforms of p and v are given by 

t 



Rfi(z) = z and R Dt{u) {z) 



and accordingly 



1 - tz 

t 



Therefore, the Cauchy transform G fl SD t (i'){ z ) * s the root of the cubic equation 

(13) tg 3 ~g 2 (l + tz) + gz-l = 0. 

First, in order to show uniqueness, we will show that Equation (|13p has only one 
root for which lim^i^^ zg(z) = 1. This follows from the observation that if g%, gi 
and 33 are the roots of Equation (fT3")) then gigigz — \ and g\ + 32 + 33 = \ + z. 
Combining both identities results in 

(14) g\gi + gxgt + j = -9152 + zg x g 2 

Thus if two of the roots, say g\ and 32, are such that Umui^^ zg\(z) = 1 and 
limui^oo zgi{z) = 1, then Equation (11411 would lead to the contradiction that -? = 0. 
It is known, see [3j Corollary 2, Corollary 4, and Proposition 5], that the free 
convolution of a compactly supported measure with a semicircle law has a smooth 
and bounded density. Thus by picking the right root and then using the Stieltjes 
inversion formula (|12p we get 

-- limS(G / _ lffl _ Dt ( y) (2; + iy)) = f(x;t) 

7T ylO 

where f(x;t) is given by Equation pip . 

Second, since we know in advance that the free additive convolution of two 
probability measures with compact supports in K has a compact support in R (c/. 
[10j). then we can find the support of p, ffl DtW) by identifying when f{x;t) = 0. 
This last identity leads to Equation (jlOp . 

The left hand side in Equation (|10p is positive at x = and Equation (fTOj) 
has two real roots and two complex conjugate roots. We prove the existence and 
uniqueness of its real roots si(t) and S2(t) as follows. Substitute x = y + + t in 
Equation (jTTJJ) . Hence, we get another quartic equation in y which reads as 

< 15 » - (f - 15? + 3 '*) + (l - 8 '0 » - + 6 (2 ) " 3 + " 4 ^ » 
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By Descartes— Euler theorem PQ, Equation (fl"5|) has two real roots and two com- 
plex conjugate roots (and correspondingly Equation (flO|) docs) if and only if the 
cubic equation in z that is given by 

(16) - (j-8* 3 ) 2 + (80 + 48i 4 ) z - (i + 12i 2 ) z 2 + z 3 = 

has one nonnegative real root and two complex conjugate roots. To show this 
we substitute z = u + ( jp- + 4< 2 ) in Equation (JTSJ) which in its turn reduces to 
another cubic equation. By Cardano's theorem, the discriminant of the new reduced 

_ 64 (l+27t 4 ) 3 . , . . . . . 

equation D — — 277 s — 1S nonn egative and so it has one real root and two complex 
conjugate roots and so does Equation (| 16[> . 

Now, it is left to show that this real root of Equation (| 16[) is nonnegative. This 
is true since the product of the three roots of Equation (fi75)) is equal to (| — 8t 3 ) 
which is nonnegative and consequently the real root is also nonnegative. 

We can see from Equation (fT0|) that si(t) and S2(t) are continuous and hence 
uniformly bounded on any compact subset of M. This completes the proof. □ 

In Figure 1, we show the graphs of f(x; t) for t — 1,2. 




Remark 3. If t = 0, then f(x; 0) = 7-5: %/4 ~ % 2 l[-2,2]( a; ) since fj. ffl Sq = fx. In case 
t < 0, since fi ffl D t (v) = ffl D- t (v)), then it follows that fi ffl D t (v)(dx) = 

f(—x; —t)dx. 

We will also need the following. 

Theorem 2. [51 Proposition 4.3.9] Let {A„} and {B„} be two sequences of Hermit- 
ian random matrices and {U n } be a sequence of random matrices with the uniform 
distribution on the unitary group U(n) . Suppose that U„ is independent o/(A„,B„) 
for all n > 1 . If there exist two compactly supported probability distributions \x and 
v such that 

/ia„ — ► H asn-KX) a.s. ana /ie„ — > v as n — + 00 a.s., 

then 

v ra 

MA„+u*B n u n — y ft ™ v as n — > 00 a.s. 
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In particular, if A„, B„ are independent, and the distribution of B„ is unitarily 
invariant for all n, i.e., if B„ and U„B„U* have the same distribution for all 
unitary matrices U„, then Theorem [2] implies 

v m 

MAn+Bn — ► [imv as n — > oo a.s. 

Since the distributions of G(n) and Wishart(n,p n ) matrices are unitarily invariant, 
we get the following. 

Corollary 1. 

(i) For n > 1, let A n be Wigner(n) and B„ be G(n) suc/i i/iai A„ and B„ 
are independent. Then for all f 6 1, 

v 

A 1 A„+tB„ — ■* 7o,i+t 2 as ra —> oo a.s. 

(ii) For n > 1, let A„ fee Wigner(n) and B„ 6e Wishart(n,p„) such that A„ 
and B„ are independent. 7/lim n _ +00 2^ = 1 7 f/ien /or a/Z f £ K, 

v 

/j,A„+tB„ — * as n —> oo a.s. 
where [it has the probability density function given by Equation ill]) . 

Now, Proposition Q] follows directly from Corollary [1] part (i). Proposition 
follows easily from Corollary Q] part (ii) since f(x;t)dx has a bounded support that 
is uniformly bounded in t £ [—2, 2] as shown in Proposition [3J 

3.3. Comments on limits with respect to one index. Here we remark about 
additional information about limits with respect to one of the indexes that can be 
extracted from the proof of Theorem [TJ 

3.3.1. Formula ([9]) identifies the limiting spectral distribution of finite dimensional 
random block-matrices as the size of the blocks goes to infinity. For instance, 
consider the sequence of k x k random block-matrices {S n }, for a fixed k G N, 
in which the diagonal blocks are made of A ra 's and all the other blocks are made 
of B n 's. This random block matrix is studied in [S] among other things, using 
algebraic manipulations of block matrices. Hence, for each n we can write S n as 
I (g) A„ + W gj B„, where I is the k x k identity matrix and W is the k x k 
non-random matrix whose entries are O's on the diagonal and one's elsewhere. By 
induction on k, one can easily find that /iw = ^jt^-i + \$k-i- If {A„} and {B„} 
satisfy condition (JTT]) in Theorem [TJ then it follows from formula ^ that 

v k — 1 1 
jU§ — ► — ; — 1, .) + -ripik — 1, .) as n — » oo a.s. 
k k 

3.3.2. The following result is a generalization of [7j Theorem 4] , where it is assumed 
that the entries of take values — 1 and 1 with probabilities \ , and A and B are 
non-random matrices. 

Proposition 4. (The SS-Law) Let {Wfe} be a sequence of random matrices such 
that Wfe is Wigner(fc) for all k. Let also A and B be two n x n Hermitian com- 
muting random matrices which have eigenvalues {ai,--- ,a„} and {Pi,-- - ,P n }, 
respectively. Then 

1 " 

(17) /j B „ fc — > - y"7a,,/3? asfc->oo a.s. 

3=1 
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Proof. Since A and B are Hermitian and commute, then there is a unitary matrix 
U such that 

A = U Diag(ai , . . . , a n ) U* and B = U Diag(/?i ,...,/?„) IT . 
It follows directly from equation ([8]) that 

limjfe^ootrnkCB^) = / R £E?=i(« J -+t&r7o 1 i(dO 

= iE;=i/ R * m 7 Q ,,^(d*) 

where the last equality follows by changing of variables. 

Since the probability measure on the right-hand side of (|17p is a finite mixture 
of semicircle laws, then it has a compact support in HL Therefore, convergence of 
moments implies weak convergence. □ 
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